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We investigate theoretically electron spin states in one dimensional (ID) and two dimensional (2D) 
hard-wall mesoscopic rings in the presence of both the Rashba spin-orbit interaction (RSOI) and the 
Dresselhaus spin-orbit interaction (DSOI) in a perpendicular magnetic field. The Hamiltonian of 
the RSOI alone is mathematically equivalent to that of the DSOI alone using an SU(2) spin rotation 
transformation. Our theoretical results show that the interplay between the RSOI and DSOI results 
in an effective periodic potential, which consequently leads to gaps in the energy spectrum. This 
periodic potential also weakens and smoothens the oscillations of the persistent charge current (CC) 
and spin current (SC) and results in the localization of electrons. For a 2D ring with a finite width, 
higher radial modes destroy the periodic oscillations of persistent currents. 
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I. INTRODUCTION 

In recent years, the spin-orbit interaction (SOI) in low- 
dimensional semiconductor structures has attracted con- 
siderable attention due to its potential applications in 
spintronic devicesi^ There are two types of SOI in con- 
ventional semiconductors. One is the Dresselhaus spin- 
orbit interaction (DSOI) induced by bulk inversion asym- 
metry^ and the other is the Rashba spin-orbit inter- 
action (RSOI) induced by structure inversion asymme- 
try4£ The strength of the RSOI can be tuned by external 
gate voltages or asymmetric doping. Recently, the intrin- 
sic spin Hall effect (SHE) in a spin-orbit coupled three- 
dimensional p-doped semiconductor- and in a Rashba 
spin-orbit coupled two-dimensional electron system^ was 
predicted theoretically. It provides us a possibility to 
generate the spin current (SC) electrically without the 
use of ferromagnetic metal or a magnetic field. 

Recently advanced growth techniques make it possi- 
ble to fabricate high quality semiconductor rings£ A 
quantum ring exhibits the intriguing spin interference 
phenomenon because of its unique topology. The per- 
sistent CC in mesoscopic rings threaded by a magnetic 
U(l) flux has been studied extensively, neglecting the 
spin degree of freedom of the electronAiSiii It has been 
experimentally observed both in a gold ring 1 - and in a 
GaAs-AlGaAs ring 13 - using standard SQUID magnetome- 
try. As for the persistent SC, the SC-induced electric field 
that was predicted by several authors^i^ may contribute 
to the successful measurement of the persistent SC in 
mesoscopic rings in future. The electronic structures and 
magneto-transport properties of ID rings with the RSOI 
alone have attracted considerable interest .ILiLiS Since 
the strength of the DSOI in thin quantum wells is com- 
parable with that of the RSOIfiS one should consider 
both of the SOI's in low dimensional structures. A few 
works have been done on the effects of the competition 
between these two types of SOI on the transport proper- 
ties of 2DEGi 20 i 21 i 22 The effects of the interplay between 
the RSOI and DSOI on the spin states and persistent 



currents (CC and SC) in mesoscopic rings are highly de- 
sirable. 

In this paper, we investigate theoretically the spin 
states and persistent CC and SC in mesoscopic rings un- 
der a uniform perpendicular magnetic field in the pres- 
ence of both RSOI and DSOI. We find that the persis- 
tent CC and SC, charge density distribution, and local 
spin orientation are very sensitive to the strength of the 
RSOI and DSOI. The interplay between the RSOI and 
DSOI leads to an effective periodic potential. This po- 
tential has significant effects on the physical properties of 
mesoscopic rings, e.g., the energy gaps, the localization of 
electrons, and weakening and smoothening of persistent 
CC and SC. Five different cases are considered: (1) a ID 
ring with RSOI alone; (2) a ID ring with DSOI alone; (3) 
a ID ring with RSOI and DSOI of equal strengths; (4) a 
ID ring with RSOI and DSOI of different strengths; (5) 
finite-width effects on the energy spectrum, charge den- 
sity distribution, the persistent CC and SC. The eigen- 
states of cases one and two are analytically solved and can 
be connected by a unitary transformation. The paper is 
organized as follows. In Sec. |H] the theoretical model 
is presented. The numerical results and discussions are 
given in Sec. IIIII Finally, we give a brief conclusion in 
Sec. El 



II. THEORETICAL MODEL 

A. Hamiltonian 

In the presence of both RSOI and DSOI, the single- 
particle Hamiltonian for an electron in a finite- width ring 
(see Fig. ^b)) under a uniform perpendicular magnetic 
field reads 

h 2 k 2 

H = ^— - + a(a x k y - a y k x ) + [3{a x k x - a y k y ) 
+ lg*^ B Ba z +V(r), (1) 
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leads to a non-hermitian Hamiltonian in the presence of 
RSOI or DSOI— We can obtain a hermitian Hamilto- 
nian including both RSOI and DSOI following Ref. I2M 
The stationary Schrodinger equation is 



iJ* = E^. 



(6) 



The dimensionless ID Hamiltonian including both RSOI 
and DSOI reads 



FIG. 1: Schematic diagrams for ID ideal ring (a) and 2D 
hard- wall ring (b). 



where k = —iV + eA/h. A(r) — B/2(—y,x,0) is the 
vector potential, m* is the electron effective mass. The 
fourth term describes the Zeeman splitting with Bohr 
magneton /j-b = eh/2m e and the effective g factor g*. 
<Ji(i = x,y,z) are the Pauli matrices, a and (3 specify 
the RSOI and DSOI strengths, respectively. V(r) is the 
radial confining potential, 



V(r) = 



7*1 ^ r ^ T2 

00 otherwise 



(2) 



where n and are the inner and outer radii of the ring, 
respectively. 

In the following we take the average radius a — (r± + 
r2)/2 as the length unit and Eq — ti 2 /2m*a 2 as the en- 
ergy unit. The dimensionless Hamiltonian in the polar 
coordinates becomes 



H = 



H k + V(r) + gb/2 /3k + + iak_ 
(3k_ - iak+ H k + V(r) - gb/2 



(3) 



where H k = (e r k r + e v k v ) 2 is the dimensionless kinetic 
term; k± = k x ± ik y = e ±l(p (k r ± ik v ), with k r = —i-j^p 
and k v — + fr; b = heB/m*Eo is the dimension- 

less magnetic field; a(/3) = a(/3)/E a specifies the di- 
mensionless RSOI (DSOI) strength; and g = g*m* /2m e 
is the dimensionless g factor. 

The wavefunction ^(r) of an electron in the ring can 
be expanded as 



(4) 



H 



i— 

dip 

+ /3 2 
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1 O r - 
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sin zip - 

2 Y 



<r<p(-<P) 



(7) 



cos tpoy — sin ipa x 



where ay — cos ipa x + sin ipa y , cr ¥ 

<I> = Bira 2 is the magnetic flux threading the ring, and 
<I>o = h/e is the flux unit. Notice that there is a peri- 
odic potential term sin 2(f) induced by the interplay 
between the RSOI and DSOI. 

We introduce a vector function S(r) to describe the 
local spin orientation of a specific eigenstate in a ID 
ring: 



S(r) = v&tsVfr = ^s x ^e x 



¥s z *e z . (8) 



When the coupling strength a or (3 vanishes, the sin 2(p 
potential accordingly disappears and the analytical solu- 
tion to the eigenvalue problem is available (see Appendix 
for details). Generally we have to solve Eq. iJBJ numeri- 
cally when a/0 and j3 7^ 0. 



B. Persistent currents 

The charge density operator and the charge current 
density operator are 



p(r') — —eS(r' — r) 
J'c(r') = ^[p{r')v + vp(r'))}, 



(9) 



where r' refers to the field coordinates and r the coor- 
dinates of the electron. We can also introduce the spin 
density operator and spin current density operator— as 



R n (r)Q m {ip)xa(s z ) 



■ sin 



717T 



[r-ri) 



^■Xcr{s z ), (5) 



where d = r-i — 1\ is the width of the ring and Xa{, a z){p = 
±1) are the eigenvectors of s z . 

Most previous theoretical studies on mesoscopic rings 
are based on the Hamiltonian of a ID ring (see Fig.^a)), 
which can be obtained by simply disregarding all the 
terms proportional to derivatives with respect to r in 
the 2D Hamiltonian. But this conventional procedure 



S(r') = -&6(r'-r) 
j s (r') = ~[S(r')v + vS(r% 



(10) 



where & = a x e x + cr y e y + a z e z is the vector of the Pauli 
matrices. The charge current density and spin current 
density can be obtained by calculating the expectation 
values of the corresponding operators: 



j c (r / ) = <*|jc|*> = -eRe{*+(r , )*'*(r / )} 
j s (r') = m'j s I*) = Re{&(r')v'm{r')} , 



(11) 
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TABLE I: Parameters used in our calculation are from Ref.126 





m*(m e ) 


9' 


GaAs 


0.067 


-0.44 


InSb 


0.014 


-51 



where vE'(r) is the wavefunction of an electron in the ring. 
For convenience, we note r', v' as r, v hereafter. 

The (^-component of the velocity operator associated 
with the Hamiltonian in Eq. JQl is 



d eBr a 



(12) 



im*r dtp 2m* h 
The azimuthal (spin or charge) current can be defined 



,10 



^ p2tt rr 2 

I v = — dp J drj v (r). 



(13) 
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FIG. 2: (a) The energy spectrum for ID GaAs rings, a = 
P = 0, g = -0.01474; (b) The energy spectrum for ID InSb 
ring, a = (3 = 0, g = -0.357. In Fig. Ha) and Fig. 0» the 
solid (dashed) lines denote spin-up (spin-down) levels. 



We ignore the Coulomb interaction between electrons 
in this work. At the low temperature, N electrons will 
occupy the lowest TV levels of the energy spectrum. The 
total (charge or spin) current is the summation over all 
occupied levels^ 

For a ID ring, the eigenstates could be expanded in the 
basis set Q m (p) — exp(irri(p) / \/2tt which is much simpler 
than that in Eq. Q, and we can get the azimuthal com- 
ponent of the velocity operator in a ID ring by specifying 
the variable r as the constant a in Eq. I|12|l . 

Most of the previous investigations of the persistent 
CC in mesoscopic rings are based on the well-known for- 
mula /„ = —dE n /d$>, where I n denotes the contribution 
to the persistent CC from the nth state and $ is the 
magnetic flux through the vmgm^^^ In this paper we 
calculate the persistent CC and SC via Eq. (|TT|) . Note 
that for a ID ring our results are identical with those 
obtained from the formula /„ = —dE n /d&. 



(a) (b) 




III. RESULTS AND DISCUSSION 

We show the energy spectrum of a ID ring in Fig.[2for 
different g factors. The relevant parameters of the mate- 
rials used in our calculation are listed in Table [i] With- 
out the spin-orbit couplings, the g factor accounts for the 
spin splitting. For a material with large g factor such as 
InSb, although the parabola behavior of the energy lev- 
els as functions of the magnetic fields is still retained, the 
periodicity of the energy spectrum is severely broken by 
the g factor, especially at large magnetic fields. For a 
material with a small g factor, e.g. GaAs, the Zeeman 
splitting is quite small even in a rather large magnetic 
field. 



FIG. 3: (a) Energy spectrum of ID ring in the presence of 
RSOI alone, where the solid lines (dashed lines) denote the 
spin-up (spin-down) levels; (b) The projection of S(r)f onto 
the x-y plane; (c) Local spin orientation for all the spin-up 
levels S(r)f; (d) Local spin orientation for all the spin-down 
levels S(r)f . a = 1, = and g = 0. 



A. ID ring with RSOI alone 

As shown in the Appendix, the electron states in a 
ID ring with RSOI alone under a uniform magnetic field 
including the Zeeman splitting can be solved analytically. 

The energy spectrum in the presence of RSOI alone is 
plotted in Fig.^a). When the g factor is set to be zero, 
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the energy of each level is given by 



E n ,<T = n 



b a a 
4 + 2 ~~ 2cos< 



tan 2 9 



(14) 



where 9 = arctan(a) (Eq. (|A.11|) or Eq. I|A.12|) represents 
the corresponding eigenstate so long as 9 n ^ a is replaced 
by 9). From the above expression we find that the spin- 
up and spin-down levels with the same quantum number 
n are separated in the b axis by 4(1/ cos 9 — 1), and both 
of the spin-up and spin-down levels are pulled down by 
tan 2 9/ A compared to the case without SOI. 

For the RSOI alone case, the local spin orientation for 
all the spin-up states along the ring is described as 

S(r)f 

= — — [sin(— #)(cos ipe x + siinpe v ) + cos(—9)e z ] . (15) 
Ana 

The local spin orientation for all the spin-down states 
S(r)^ is opposite to S(r)^ as shown in Fig. |3{c) and 
Fig- Eld). In this case, the oblique angle 9 = arctan(a) 
becomes independent from the quantum number n and 
the magnetic field b. It means that the local spin ori- 
entations for all the spin-up (spin-down) states are the 
same. When there is RSOI alone, the local spin orien- 
tation exhibits rotational symmetry for either spin-up or 
spin-down states. 

The contributions to the persistent CC and SC from 
each level can be easily obtained explicitly: 



b a a 

n-\ 1 

4 2 2 cos 9 



b a a 

n^ 1 

4 2 2cos6» 



CTCOSf 



(16a) 
(16b) 



In,a {Ina) denotes the contribution to the persistent CC 
(SC) from the eigenstate ^ n ,<j- I n ,a {In,a) is in the units 
of 2Eo/&o (Eq/2it). We should notice that these expres- 
sions deduced from operators coincide with the formula 
I n = -dE n /d®. 

Splettstoesser et al. analyzed the persistent CC in- 
duced by the magnetic flux in the ID ring with both 
RSOI and a impurity potential*^ They demonstrated 
that the strength of the RSOI can be extracted from the 
dependence of the persistent CC on the magnetic flux. 
The number of electrons N in their work was assumed 
to be large enough. We focus on the case in which there 
are few electrons in the ring (see Fig. . We find that 
the persistent CC is a periodic function of b exhibiting 
many linear segments with a slope ratio of —1/4 which 
can be easily deduced from Eq. H16al) . The periodicity of 
the persistent CC for an arbitrary N is 4, the same as 
that of the energy spectrum. For an arbitrary even num- 
ber of electrons N = 2m, the jumping amplitude is 1/2 
and the neighboring two jumps are shifted along the CC 
axis. For an odd number of electrons N = 2m + 1, there 




FIG. 4: The persistent CC in a ID ring with different numbers 
of electrons N vs magnetic field b while a = 1, /3 = 0, g = 0. 
The persistent CC is in units of 2NE /& . 



are two jumping amplitudes which appear alternately at 
those points where b = 0, 2, 4, 6, . . . One jumping ampli- 
tude is (m — l/cos£? + 2)/(2m + 1), while the other is 
(m + 1/ cos# — l)/(2m + 1). When n approaches infinity, 
the above two amplitudes tend to 1/2. 

The dependence of the persistent SC on the magnetic 
field b shows interesting behavior. Turning (jumping) 
points in the persistent SC oscillation for an odd (even) 
number of electrons are caused by the crossing of lev- 
els with opposite (same) spins. When the magnetic field 
b sweeps, the persistent SC oscillation for an odd (even) 
number of electrons exhibits saw-tooth (square) wave be- 
havior and the oscillation amplitude of the persistent SC 
for an even number of electrons is bigger than that for 
a neighboring odd number of electrons especially when 
the number of electrons increases. For an odd number of 
electrons N = 2m + 1, the slope ratios of the persistent 
SC are ±cos6>/4(2m + l) alternately. For an even number 
of electrons N = 2m, the jumping amplitude is cos 9/2 for 
all jumping points. 



B. ID ring with DSOI alone 

Now we consider a ID ring with DSOI alone. The 
eigenstates of the Hamiltonian of a ID ring with a = a, 
f3 = b, and g — c can be connected to those with a — b, 
(3 = a, and g — —c by a unitary operator T (see the 
Appendix). As we will show in the next subsection, T 
represents a rotational transformation in spin space, b 
and c are set to be zero and the relationship between the 
RSOI alone case and DSOI alone case is specified as 



exp[in/i]T^ tl , 



*£ a =expH7r/4]Tt*£ T . 



(17a) 
(17b) 
(17c) 
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FIG. 5: Same as Fig. 0] but for the persistent SC. The per- 
sistent SC is in units of NEo/2n. 




(b) 



(d) 




The energy spectrum while a = 0, j3 = 1, and 5 = is 
plotted in Fig. Eta), which is exactly the same as that of 
a ID ring with RSOI alone (see Fig. Eta)), but the spin 
orientations of the corresponding eigenstates are differ- 
ent. This feature arises from the sign reversal behavior 
of a z under the unitary transformation T. 

Although the eigenstates for the RSOI alone case and 
those for the DSOI alone case can be connected by the 
unitary transformation T, the behaviors of S(r)^ and 
S(rW, are very different. In the case with DSOI alone, 
spin-up (spin-down) states with different angular quan- 
tum number n or under different magnetic field b share 
the same local spin orientation S(r)P (S(r)?). The an- 
gle between S(r)P (S(r)P) and the z-axis is —9 (it — 9). 
The local spin orientations S(r)P can be obtained by 
interchanging the x and y components of S(r)?] (see 
Eqs. Q5) and |T5|)). 



S(r)f 



47ra 



[sin(— (9)(sin ^e^ + cosv5e y ) + cos(— 6)e z ] . (18) 



This interesting feature comes from the behavior of <J x (y) 
under the unitary transformation T, i.e., Ta x ty\T* = 
T^a x ( y )T — —a y ( x y The projections of S(r)^ and S{r)P 
onto the x-y plane are very different (see Fig. Etb) and 
Fig. Etb)). F° r the RSOI alone case, the vector always 
points along the radial direction and the locus of the ar- 
rowhead is a circle (see Fig. Etb))- For the DSOI alone 



FIG. 6: (a) Energy spectrum of ID ring in the presence of 
DSOI alone, where the solid lines (dashed lines) denote the 
spin-up (spin-down) levels; (b) The projection of S(r)® onto 
the x-y plane; (c) Local spin orientation for all the spin-up 
levels S(r)f ; (d) Local spin orientation for all the spin-down 
levels S(r)f . a = 0, = 1 and g = 0. 



case, the vector varies along the ring and the cylindrical 
symmetry is broken (see Fig. Etb)). 

In the current case, the persistent CC oscillations ex- 
hibit the same behavior as those of the RSOI alone case 
(see Fig.Q} because the corresponding levels are identical 
as functions of the magnetic field b. We know that the 
contribution of each level to the persistent SC is related 
not only to the magnetic field dependence of the eigenen- 
ergy but also to the spin orientation of the eigenstate. 
Since the spin-up and spin-down levels are interchanged 
compared to those of the RSOI alone case, the persistent 
SC in the current case can be obtained by changing the 
sign of the persistent SC for the RSOI alone case. 



C. ID ring with equal strength RSOI and DSOI 

The electron energy spectra for a = /3 = l, g = 
and a = [3 = 3, g — are plotted in Figs. Eta) an d 
Etc), respectively. We can see that the energy spectra 
are spin degenerate. The spin degeneracy comes from the 
symmetry of the Hamiltonian when a — [3. In this paper, 
we use the unitary operator T (see Appendix) to describe 
the symmetry of the ID Hamiltonian. When a — (3 7^ 
and 3 = 0, THT^ = T^HT = H. If we have = 
E^, in which 'J is an eigenstate for the eigenenergy E, 
then the states T\& and T'^f are also eigenstates and 
are equivalent to each other. Thus the energy levels are 
twofold degenerate. The operator for a tp rotation around 
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FIG. 7: (a) Energy spectrum of ID ring while a = j3 = 1, 
g = 0; (b) The charge density (CD) distribution of the lowest 
single electron state in a ID ring while a = f3 = 1, g = 0, 
6 = 0; (c) Energy spectrum of ID ring while a = (3 = 3, 
<? = 0; (d) The charge density (CD) distribution of the lowest 
single electron state in a ID ring while a = f3 = 3, p = 0, 
6 = 0. The charge density is in units of e/2ira. 



the unit vector n in the Hilbert space reads 

where i and s denote the orbital and spin angular mo- 
mentum operators, respectively. It is easy to demon- 
strate that the unitary operator T can be written as 
exp[-Mrni • s/h] with m = (1/^2,-1/^2,0). Then T 
is actually a rotation operator in the spin space. For a 
quantum ring in the x-y plane, the orbital angular mo- 
mentum vector L = r x p points along the z-axis, and 
therefore ni -L = 0. Thus the unitary operator T is even- 
tually a rotation operator in the whole Hilbert space: 



T = D(m,ir) 



(20) 



That means the unitary transformation T (T<) is actu- 
ally a rotation by it (—it) around ni. Similarly, there are 
also symmetric operations corresponding to it and — it ro- 
tations around ri2 = (l/y/2, 0) for the Hamiltonian 
with a = —(3. We need to stress that these symmetric 
operations also exist for a two-dimensional electron gas 
with equal strength RSOI and DSOI (a = ±/3) and g — 0. 

Besides the spin degeneracy, it is interesting to find 
that gaps appear in the energy spectra. In order to un- 
derstand this feature, we transform the original Hamil- 
tonian to a simple form via a unitary transformation A 



as follows: 
A* 

exp[-iaf(<p)] exp[iaf(tp)] 

exp[— itt /A] exp[— iaf(<p)] — exp[— i7r/4] exp[iaf(tp)] 

(21) 

with f(<p) = sin((/3 + 7r/4). The original electron Hamil- 
tonian with a = j3 and 5 = reads 



d b a a , . 
H= - l - + l + -a r --a v (- i p) 



~Y + Y s[n2ip - (22) 

After the transformation, the Hamiltonian becomes 



H' = A^HA 



dip 



a 
Y 



■sm2(p. (23) 



That means the Hamiltonian of a ID ring with equal 
strength RSOI and DSOI and zero g factor is equivalent 
to that of a ID ring with a periodic potential alone (see 
the last term in Eq. I|23(l '). The potential height is propor- 
tional to the square of the SOI strength, and the average 
of the potential shifts down by about a 2 /2. The eigenvec- 
tors of Eq. []2.'iH are actually the periodic solutions of the 
Mathieu equation^! The energy gaps, which are propor- 
tional to the potential height, decrease with decreasing 
SOI strengths, especially for higher gaps (see Fig. da)). 
When the strengths of SOI are fixed, the higher energy 
gaps are narrower than the lower ones because there is 
less influence from the potential. 

When only one type of SOI (RSOI or DSOI) exists, 
the charge density distribution will be constant along 
the ring. But the charge density distribution becomes 
localized along the ring when both RSOI and DSOI are 
taken into account. This localization arises from the ef- 
fective periodic potential, whose height is determined by 
the product of the strengths of RSOI and DSOI (see 
Eq. Therefore large SOI strengths lead to strong 

electron localization. The absolute value of the charge 
density exhibits maxima at the valleys of the su^^ po- 
tential (ip — Ztt/A or tp — 7 'it /4) and minima at the 
peaks of the sin2<^ potential (<p = 7r/4 or <p = 57r/4) 
(see Fig. □» and Fig. W[d)). 

In the spin degenerate case, we cannot define the lo- 
cal spin orientation S(r) for an eigenenergy level and the 
persistent SC for an odd number of electrons because of 
the uncertainty of the eigenvectors. The persistent CC 
oscillation for an odd number of electrons N = 2m + 1 
is simply the arithmetic average of that for N = 2m and 
that for N = 2m + 2 in a spin degenerate case. Thus we 
compare the persistent current (CC and SC) oscillations 
in the three degenerate cases: a — (3 = 0, a — (3=1, 
and a = (3 = 3 only for even numbers of electrons (see 
Fig. EJ. The persistent SC is zero in the two degenerate 
cases for even numbers of electrons. The sin 2ip potential 
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FIG. 8: The persistent CC (solid lines) and SC (dashed lines) 
with different even numbers of electrons TV vs magnetic field 
in the degenerate cases (a) a = j3 = g = 0; (b) a = /3 = 1, 
g = 0; (c) 5 = P = 3, g = 0. The persistent CC (SC) is in 
units of 2NE /$ (NE /2w). 



in Eq. Q accounts for the natter magnetic dispersion as 
well as gaps in the energy spectrum when a = f3 =/= 0. 
Since the contribution to the persistent CC from an en- 
ergy level is actually determined by the dependence of 
the energy level on magnetic field, the oscillation of the 
persistent CC for the case a = (3 ^ is smoother and 
smaller than that for the case a = f3 = 0. The interplay 
between the RSOI and DSOI smoothens and weakens 
the persistent CC oscillation most obviously when the 
Fermi level locates near the lowest gap (see the panels 
labeled N = 4 in Fig. \^b) and Fig. Q^c)). While the 
number of electrons increases, the oscillation of the per- 
sistent CC becomes sharp again since the higher gaps 
become smaller. This smoothening and weakening effect 
can even be found again for a large number of electrons 
when the SOI strengths increase (see Fig.JSfc)). 



D. ID ring with different strength RSOI and DSOI 







FIG. 9: (a) Energy spectrum of ID ring while a = 2, /3 = 1, 
g — 0; (b) The charge density (CD) distribution of the lowest 
single electron state in a ID ring while 5 = 2, P = 1, g = 0, 
6 = 0; (c) Energy spectrum of ID ring while a = 4, j3 — 3, 
g — 0; (d) The charge density (CD) distribution of the lowest 
single electron state in a ID ring while a = 4, f3 — 3, g — 0, 
6 = 0. The charge density is in units of e/2na. 
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Generally, the symmetry of the Hamiltonian shown in 
the previous subsection no longer exists when \a\ ^ |/3|, 
even for g = 0. We show the electron spectra for a = 2, 
j3 = 1, g = and a = 4, ft = 3, g = in Fig. E|a) and 
Fig- He), respectively. The energy gaps increase with in- 
creasing SOI strengths. But the spin splitting in the two 
spectra is quite different. It is interesting to note that the 
energy spectrum becomes spin degenerate again when the 
two SOI strengths are tuned to proper values even though 
they are different. Figs. |^b) and^d) show the distri- 
bution of charge density for different strength RSOI and 
DSOI. The electron is localized along the ring due to the 
periodic potential ^ sin 2ip, The electron density distri- 
bution becomes more localized with increased potential 



FIG. 10: The persistent CC in a ID ring with different num- 
bers of electrons N vs magnetic field b while a = 2, f) = 1, 
g = 0. The persistent CC is in units of 2NE /$o- 



height, i.e., the product of the strengths of RSOI and 
DSOI (see Eq. Q). 

The persistent CC and SC are plotted in Fig. ^| and 
Fig. El respectively. We find that the oscillations of 
the persistent CC and SC become smooth and weak due 
to the gaps in the energy spectrum, especially when the 
Fermi level locates near the largest energy gap (N = 4). 
The persistent current (CC or SC) oscillation no longer 
consists of linear segments since the parabolic behavior of 
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b b 



FIG. 11: Same as Fig. QJ3 but for the persistent SC. The 
persistent SC is in units of NEo/2ir. 




the energy levels disappears due to the periodic potential 
^ sin2<£> in the ID Hamiltonian (see Eq. Q). 

The local spin orientation S(r) also reveals the inter- 
play between the RSOI and DSOI since S(r) R is quite 
different from S(r) D . According to Eq. J7J), the interplay 
between the RSOI and DSOI is divided into two parts, 
i.e., w<r r — v) m the kinetic term and the periodic 

potential ^ sin %p. The first part makes the direction of 
the local spin orientation vary along the ring, and the sec- 
ond part leads to the electron localization (see Figs. ll2f a) 
and I12f b11. When the SOI strengths increase, the spin 
orientation exhibits rapid variation due to the enhance- 
ment of the interplay between the RSOI and DSOI. 

In Fig. ll2f cL we plot the persistent SC as a function of 
RSOI strength a and DSOI strength (3 for a fixed number 
of electrons N — 8 and magnetic field 6 = 2. The contour 
plot shows interesting symmetry. It is symmetric (anti- 
symmetric) with respect to the lines a = and j3 = 
(a = ±/3). From this figure we find that the maxima 
and minima of the persistent SC occur while only one 
of the two types of the SOI exists. That is because the 
effects of the RSOI and DSOI on spin splitting tend to 
cancel each other. The persistent SC becomes zero when 
the strengths of the RSOI and DSOI are equal to each 
other (a = ±/3). This corresponds to the spin degen- 
erate case discussed before. Besides the two orthogonal 
lines (a = ±/3) on the a-j3 plane there are many circle- 
like closed curves on which the persistent SC disappears. 
These zero-SC lines intersect the a axis at those points 
(±Vm 2 — 1, 0) and the (3 axis at (0, ±y/m 2 — 1) where 
m = 1,2,3,.... The persistent SC disappears because 
the energy spectrum becomes degenerate again when the 
strengths of RSOI and DSOI are tuned to proper val- 
ues even though they are not equal. The contour plot 
of the oscillation amplitude of the persistent CC for a 
fixed number of electrons N = 8 is shown in Fig. I12f d^l . 
The oscillation amplitude as a function of a and (3 is 



FIG. 12: (color online) (a) Local spin orientation S(r) for the 
lowest spin-up level while a = 2, f3 = 1, g = 0, 6 = 2; (b) 
Local spin orientation S(r) for the lowest spin-up level while 
a = 4, (3 = 3, g = 0, b = 2; (c) The persistent SC in a ID ring 
with different RSOI and DSOI strengths when the magnetic 
field b is 2; (d) The oscillation amplitude of the persistent 
CC in a ID ring with different RSOI and DSOI strengths. 
In Fig. IT3T c1 and Fig. THtd) we set g = and N = 8. The 
persistent CC (SC) is in units of 2NE /$ (NE /2tt). 

symmetric with respect to the lines a = and (3 = 
and a = ±/3. When only one type of SOI appears, the 
maximum of the persistent CC oscillates with increased 
SOI strength. When both types of SOI are included, the 
maximum of the persistent CC decays since the interplay 
between the RSOI and DSOI leads to a periodic poten- 
tial along the ring which results in the gaps in the energy 
spectrum, consequently smoothening and weakening the 
oscillation of the persistent CC. 



E. Finite width effects 

Now we turn to consider a mesoscopic ring with finite 
width. When a 2D ring is thin enough, its character- 
istics are almost the same as those of a ID ring since 
the second radial levels are too high to be occupied and 
the compressing effect of the magnetic field on the radial 
wave function is negligible. Here we consider the wide 
ring case. 

The energy spectrum for a 2D ring with a finite width 
while a = 2, fj = 1 and g = is plotted in Fig. H3f aL 
The second radial mode can be seen in the top of the en- 
ergy spectrum. The compressing effect of the magnetic 
field on the radial wave function accounts for the increase 
of the energy with increasing magnetic fields. Like the 
ID case, the interplay between the RSOI and DSOI leads 
to an effective periodic potential (see Eq. (|A.2|) L result- 
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(a) 



(b) 





FIG. 13: (a) Energy spectrum for 2D hard wall ring with 
width d = 1.33 while a = 2, /? = 1, g = 0; (b) The charge 
density (CD) distribution of the lowest single electron state 
in a 2D hard wall ring with d — 1.33 while a = 2, j3 = 1, 
g = 0, 6 = 0. The charge density is in units of e/2na 2 . 




FIG. 14: (a) The perisistent CC in a 2D ring vs magnetic 
field b while a. = 2, /3 = 1, g = 0, the width d = 1.33 and 
N — 16; (b) Same as (a) but for the persistent SC; (c) Same 
as (a) except N = 20; (d) Same as (b) except N = 20. The 
persistent CC (SC) is in units of 2NE /3> {NE /2tt). 



IV. SUMMARY 



We have conducted a theoretical investigation of the 
spin states and persistent CC and SC in mesoscopic rings 
with spin-orbit interactions. We have demonstrated theo- 
retically that the Hamiltonian of the RSOI alone is math- 
ematically equivalent to that of the DSOI alone by a uni- 
tary transformation T. This property results in the de- 
generate energy spectrum for equal strength RSOI and 
DSOI. The interplay between the RSOI and DSOI leads 
to an effective periodic potential ^ sin 2tp. This peri- 
odic potential results in gaps in the energy spectrum, 
and smoothens and weakens the oscillations of the per- 
sistent CC and SC. The charge density and the local spin 
orientation S(r) are localized along the ring due to the 
effect of the periodic potential. Higher radical modes be- 
come involved as the ring width increases, destroying the 
periodicity of the persistent CC and SC oscillations. 
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ing in energy gaps that depend on the magnetic field. 
The charge density distribution of a single electron in a 
2D ring at b = is shown in Fig. H3T b). We find that 
the electron probability is localized due to the sin 2<p po- 
tential arising from the interplay between the RSOI and 
DSOI. 

The persistent CC and SC in the 2D wide ring as a 
function of magnetic field are plotted in Fig. [21 The 
number of the electrons is tuned to detect the effect of the 
second radial mode. When we consider the contribution 
of the lowest 16 levels, the oscillations of the persistent 
CC and SC have a profile similar to that in the ID ring 
because the second radial mode has not been involved 
yet. But when we increase the number of the electrons 
to 20, the quasi-periodicity of the persistent CC and SC 
as functions of the magnetic field b is destroyed. This 
feature arises from the contribution of the second radial 
mode. 



APPENDIX: THE HAMILTONIAN AND 
AVAILABLE ANALYTICAL SOLUTIONS 



The dimensionless Hamiltonian for a 2D ring reads 



H 2D =H k +H R + H D + H z + V(r), (A.l) 



where the kinetic term = [e r k r + e lfi k ip ) 2 , the Rashba 
term Hr = a{a r k v — <r v k r ), the Dresselhaus term Hp = 
(3[a r (—<p)k r — Gip{—<p)km\, and the Zeeman term Hz = 
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gba z /2. V(r) is the radial confining potential. 



H, 



2D 



dr 



ld_ 

r dr 



V(r) 



+ [ao r — j3a v {—ip)\ 

-<i 

-acr ip + /3oy(— ip)] 



i - 1 



. _ 
lr 2r 
d 2 Id 
dr 2 r dr 



+ V(r) 
-aa v + (3a y 



k r - 
2r 



a 

2 ar 



+ 2 9b(7z 



dr z r dr 



k^ 

ken ~t~ 



[~aa v + /3a r (-<p)] 




-aJ-ip) 



a 2 + 1 a/3 



1 



sin 2ip + —gba 



(A.2) 



Specifically we write H2B = Hq + H±, where Hq = 
— ■§pz — ^-§^ + V(r). The correct ID Hamiltonian H can 
be obtained by evaluating the expectation of H\ in the 
lowest radial mode of Hq 2 ^ In the limit of a very narrow 
ring, we can set r to be a constant value (r — 1) and 
the following equation will hold for an arbitrarily given 
confining potential V(r): 



d 1 

(Pal TT + ^T \Po) = 0. 
dr 2r 



(A.3) 



Here po is the lowest radial mode for V(r). Now we can 
write the ID Hamiltonian explicitly. We get 



H 



.d_ 

dip 
a 2 +/3 2 



b a 
4+2^" 

I — — sin 2ip ■ 



o cr v(- ( ^) 



(A.4) 



A unitary operator T = 



is 



exp[— in /A] 

exp[z7r/4] 
defined, and we have = T _1 = —T. By applying this 
unitary operator, the Hamiltonian becomes 



THT t 



d 



b f3 
4 + 2 ar 



dip 

d 2 +fi 2 af3 . 



1 



sin 2<p - -g6cr 2 



(A.5) 



Thus the Hamiltonian in which a = a, (3 = b, g = c is 
mathematically equivalent to that in which a = b, J3 = a, 
9 = -c 

The Hamiltonian in matrix form is 



H 



Hn H12 
H21 H22 

. d 



where 



dip 



■9b/2, 



H 



22 





. d 

dip 


ae lv (- 


. d 

dp 








+ V 



b 1 

b 1 

4 + 2 



i/3e- lv 



dip 
,d_ 

dip 



gb/2. 



b 1 

4 + 2 

b 1 

4 ~ 2 

(A.6) 



To solve the equation H'Q = E^>, we expand the wave- 

- ^ (IZ) e.(H where 
Q m (ip) = exp [im^]. The secular equation becomes 



function $ as $ 



b rn+ ia(m + I + 5) + i6 m _i/3(m + \ - \) 
= [E-(m+ |) 2 - §6/2] a m 

a m _i«(m + I - |) - ia m+1 /3(r7i + | + 5) 
= [£-(m+|) 2 +§V2] & m 



(A.7) 



Generally, we can write the Hamiltonian in an infinite 
quintuple diagonal matrix form based on Eq. l|A.7j) . 
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a_i 
b-i 
a 
bo 
ai 
bi 



gb 

2 


















0-1 






= £ 
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6 




01 




6i 



(1-1) -f 










(!) 2 + 



a( 








m+h) (i+i) 2 +f 






(!) 



o 



2 _ gb 

2 









(! + i) 2 - 



g6 



(A.8) 



We consider three different cases. For the first and sec- 
ond cases, in which the quintuple diagonal matrices are 
reducible, analytical solutions can be obtained. 



(to + 6/4) 2 + gb/2 a(m + 6/4 + 1 /2) 
a(m + 6/4+l/2) (m + b/4 + I) 2 - gb/2 



and 



= -. — [sin(7r - O n \ )(cosipe x + sin^e„) + cos(7r - O n \ )e z ] . 
47ra 



E 



b m +i 



bm+l 

(A.9) 



(A.14) 



The eigenvalues are 



Ei 



a 



b G 

n H 1 

4 2 2cos<?„, 



where tanfl n , CT = »°l/4+^2-g&/2 ■ Thc corresponding 
eigenvectors are 



tan 2 6>„ !(T 56 



b) a = 0, /? ^ 0; 

"(m + 6/4 + l) 2 + .gb/2 i/?(m + 6/4 + 1/2)' 
-i/3(m + 6/4+1/2) (to + 6/4) 2 - #6/2 



ttm+l 

6 m 



+(J 



2 COS n 

(A.10) 



= £7 



o m +i 

6m 



(A.15) 



and 



l_,(n+l/2) v f cos^e W 2 
2^ Usin^le^/ 2 , 



*(n-i/2) v ( sin%-e 



2tt 



cos %i e ^/2 



(A.11) 



(A.12) 



The eigenvalues are 
^'=( n+ i-f + 2^ 



cr 



2 tan 2 ?7„ j(T 36 



-+CT- 



2cos77„ j(T ' 
(A.16) 

where tan?7 ni(T = n +b/t-f/2+sb/2 • Tne corresponding 
eigenvectors are 



The local spin orientations for eigenstates are 

= -p— [sin(-0 n T )(cos<£e x + sin^) + cos(-fl„ !T )e z 
47ra 



1 



2tt 



and 



(A.13) 



n,l 



i(n+l/2) V 



27T 



,^l e ^/2 

-isin^e-^/ 2 



cos^e"^/ 2 



(A.17) 



(A.18) 
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The local spin orientations for the eigenstates are c) a ^ 0, (3 ^ 0. 

[sin(-r? n) T)(sin(^e :c + cos ipe y ) + cos(-rj n ^)e z ] 



Ana 



(A.19) 



and 



While both RSOI and DSOI have nonvanishing 
S( r ) n ,i strengths, we cannot reduce the infinite quintuple diago- 



U nal matrix shown in Eq. I|A.8(I into a more compact form. 

— [sin(7r - ^.^(srn^ + cos^e,) + cos(7r - Vnd )e z ] . Thug analytical so i utions do not exist. We give numerical 

(A. 20) results instead. 
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